Solitons in Bose-Einstein Condensates with time-dependent atomic scattering length 

in an expulsive parabolic and complex potential 
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We present two families of analytical solutions of the one-dimensional nonlinear Schrodinger equa- 
tion which describe the dynamics of bright and dark solitons in Bose-Einstein condensates (BECs) 
with the time-dependent interatomic interaction in an expulsive parabolic and complex potential. 
We also demonstrate that the lifetime of both a bright soliton and a dark soliton in BECs can be 
extended by reducing both the ratio of the axial oscillation frequency to radial oscillation frequency 
and the loss of atoms. It is interested that a train of bright solitons may be excited with a strong 
enough background. An experimental protocol is further designed for observing this phenomenon. 
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D I. INTRODUCTION 

-i—> 

The Bose-Einstein condensates (BECs) at nK temperature can be described by the mean field theory - nonlinear 
Schrodinger (NLS) equation with a trap potential, i.e. .the Gross-Pitaevskii (CP) equation. Recently, with the exper- 
imental observation and theoretical studies of BECs [lj], there has been intense interest in the nonlinear excitations 
of ultra-cold atoms, such as dark solitons 0, [H, 0, H, IS 0] , bright solitons [1, [H, [l(| , vortices [HI and the four- wave 
mixing [T3 |. Recent experiments have demonstrated that the variation of the effective scattering length, even including 
its sign, can be achieved by utilizing the so-called Feshbach resonance [H, [TH, EH] ■ It has been demonstrated that 
the variation of nonlinearity of the GP equation via Feshbach resonance pro vides a powerful tool for controlling the 
O . generation of bright and dark soliton trains starting from periodic waves [16[ . 

At the mean-field level, the GP equation governs the evolution of the macroscopic wave function of BECs. In the 
physically important case of the cigar-shaped BECs, it is reasonable to reduce the GP equation into a one-dimension 
■ nonlinear Schrodinger e quat ion with time-dependent atomic scattering length in an expulsive parabolic and complex 
OO ' potential [H El El 119, ML HI , 

o : 

i«-0 + 2a(t)M 2 ^-iAVV> + i7V', (!) 

\6 ' 



where the time t and coordinate x are measured in units 2/uj± and a±, a± = ^/h/muj± and ao = i/ft/mwo are linear 
oscillator lengths in the transverse and cigar-axis directions, respectively. uj±_ is the radial oscillation frequency and ojq 
is the axial oscillation frequency, m is the atomic mass, |A| = 2|wo|/w^ <C 1, a(t) is a scattering length of attractive 
, interactions (a(i) < 0) or repulsive interactions (a(t) > 0) between atoms, and 7 is a small parameter related to the 
feeding of condensate from the thermal cloud [23| . When a(t) = goexp(Xt) and 7 = 0, Liang et al. pre sent a family 
of exact solutions of ([1]) by Darboux transformation and analyze the dynamics of a bright soliton [17| ■ Kengne et 
al. investigated ([1]) with a(t) — go and 7 = A/2 and verified the dynamics of a bright soliton proposed [HI]. These 
results show that, under a safe range of parameters, the bright soliton can be compressed into very high local matter 
densities by increasing the absolute value of the atomic scattering length or feeding parameter. 

In this paper, we develop a direct method to derive two families of exact solitons of Eq. , then give some thorough 
analysis for a bright soliton, a train of bright solitons and a dark soliton. Our results show that for BEC system with 
time-dependent atomic scattering length, the lifetime of a bright or a dark soliton in BECs can keep longer times 
by reducing both the ratio of the axial oscillation frequency to radial oscillation frequency and the loss of atoms. It 
is demonstrated that a train of bright solitons in BECs may be excited with a strong enough background. We also 
propose an experimental protocol to observe this phenomenon in further experiments. 
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II. THE METHOD AND SOLITON SOLUTIONS 



We can assume the solutions of Eq. ([I]) as follows 



, , , . . , . S cosh(0 + cosfn) , „ . . a sinh(£) + 8 sinfn), 
V> - [A Q (t) + Ai(t) — ^ 44 + LBi(i) . c A ] exp(iA), 



cosh(£) + 5 cos(?7) 



cosh(£) + (5 cos(?7) 



(2) 



where 



A = ko(t) + ki(t)x + k 2 (t)x 2 , £ = Pl (t)x + qi (t), n=p 2 (t)x + q 2 (t), 



and A (t), Ai(t), B\(t), pi(t), qi(t), P2(t), Q2(t), ko(t), ki(t), k 2 {t) are real functions oft to be determined, and a,/?, 5 
are real constants. 

Substituting Eq. @ into Eq. (|TJ) , we first remove the exponential terms, then collect coefficients of 
sinh*(£) cosh J (£) sin m (?7) cos n (ri)x k (i = 0, 1,2, j = 0,1; m = 0,1,2, ■ ■ •; n = 0, 1; k = 0, 1, •••.) and separate 
real part and imaginary part for each coefficient. We derive a set of ordinary differential equations (ODEs) with 
respect to a(t), A (t), A\(t), B\{t), pi(t), qi(t), P2(t), q%(t), k (t), ki(t), k 2 (t). Finally, solving these ODEs, we can 
obtain two families of analytical solutions of Eq. (JTJ) . 

Family 1. When interaction between atoms is attractive such as 7 Li atoms, ait) < 0, the solution of Eq. (1) can 
be written as: 



where 



ipi = n[A c + A s 

A 

a(t) 

n 



(S cosh(£) + cos(?7)) 
cosh(£) + S cos(rj) 



. . a sinh(f) + 8 sinfn), 

[A ° — AA/ A exp iA + 7t - 

cosh(£) + o cos{rj) 



(3) 



k 2 (t)x 2 + k 1 fl 2 x + (2g A 2 c -kl) /ft 4 dt, 

2g oP 2A 2 c 



g^A s [pn z x-2[k 1 8+p 2 - 



(g A 2 +p 2 ) 



p 2 tt 2 x - [2p 2 fci + (pt - 9qA 



-gofl 2 exp(-27f), a 
p 2 + g (A 2 -4A 2 c ) i 



2^/goAcP2 
' 9vA 2 s +p 21 
2g Q A c A s 



9oA 2 +p 2 goAl+p 2 ' 

/A A 
-2k 2 (t)dt], k 2 (t)={±-, -tanh(Ai)}, 



and A c ,A s ,go > 0, p 2 , k\, 7 are arbitrary real constants. 

Family 2. When interaction between atoms is repulsive such as 23 Na and 87 Rb atoms, a(t) > 0, the solution of 
Eq. (1) can be written as: 



where 



V>2 = OL4 c + iA s tanh(£)] exp(zA + 7*), 



£ = ±^n 2 x + 2A s {^k 1 +g A c ) J&dt, 



k 2 (t)x 2 + k 1 n 2 x-[2g (A 2 c +A 2 s ) + k 2 } Kl 4 dt, 



(4) 



A 

a(t) = goti 2 exp(— 27*), 

and Q, k 2 (t) are the same as in (J3|), A c , A s ,go > 0, k\ and 7 are arbitrary real constants. 

The solutions ([3]) and ([4]) are new general solutions of equation {T]) which can describe the dynamics of bright 
and dark solitons in BECs with the time-dependent interatomic interaction in an expulsive parabolic and complex 
potential. In special case, it can be reduced to solutions obtained by others. For example, if k 2 (t) = — A/4 and 7 = 0, 
the solution ^ describe dynamics of a bright soliton in BECs with time-dependent atomic scattering length in an 



3 



expulsive parabolic potential, and it can reduce the solution in Ref. [T3|- If k%{t) = —A/4 and 7 = A/2, Eq. ([3]) 
describe dynamics of bright matter wave solitons in BECs in an expulsive parabolic and complex potential, and it 
can recover the solution in Ref. 18]. 

To our knowledge, the other solutions from Eqs. and have not been reported earlier. When k%(t) — ±A/4 
and 7 is a fixed value, the intensities of Eqs. ([3]) and ([4]) are either exponentially increasing or exponentially decreasing 
so the BECs phenomenon can not be stable reasonably. Thus in order to close experimental condition and compare 
our theoretical prediction with experimental results, we will only discuss and analyze Eqs. ([3]) and (|U) with k 2 (t) — 
A/4tanh(A£) and Q 2 = sech(At)/2. 

A. Dynamics of Bright Solitons in BECs 

In the following, we are interested in two cases of Eq. Q . 
(I) When a = and p 2 = 0, ipi can be written as 

/ a\A _i_ A 5 cosh (£) + cos(?7) + i/3 sin(7?) 

tpn = il[A c + A s — xexp(iA + 7t), (5 

cosh(£) + cos(^) 



where A 2 S > AA 2 c: and 



Z = ^(3A s n 2 x-2^A s k 1 f3 Jn 4 dt, 



P J A* ' 

A = k 2 {t)x 2 + k 1 n 2 x + (2g A 2 c -kf) J f2 4 di, 

_ /sech(At) . 2A C 
a(t) = — ^sech(At)exp(— 27t). 

When A s = 0, tpn reduces to the background 

ip c = A c nexp(iA + jt). (6) 
When A c = 0, ipu reduces to the bright soliton 

i/j s = A s fisech(^) exp [i9 + "ft] , (7) 

where 



i = ^A s fl 2 x~2^A s {k 1 + P2 ) J Q 4 dt, 
6 = k 2 (t)x 2 + k 1 n 2 x+[g A 2 s -k 2 } J n 4 dt. 



Thus ipu represents a bright soliton embedded in the background. At the same time, when A c <C A s satisfied 
4A 2 < A 2 and 7 small, the background is small within the existence of bright soliton. 

Considering the dynamics of the bright soliton in the background, the length 2L of the spatial background must 
be very large compared to the scale of the soliton. In the real experiment [8J], the length of the background of 
BECs can reach at least 2L = 370/im. In Fig.l, the width of the bright soliton is about 21 — 14nm [a unity of 
coordinate, Aa; = 1 in the dimensionless variables, corresponds to aj_ = (hj (mui^) 1 ^ 2 = lA/im}. So / <C L, a 
necessary condition for realizing bright soliton in experiment. From Fig. la, under the realistic experiment parameters 
in i.e., uj± — 2n x 710Hz, ujq = 2m x 70Hz. In order to cope with the experiment: the soliton move to —x 
direction, A = — 2|cjo|/^± ~ —0.197, 7 = —0.01 and 50 = 0.4nm , we can see that the lifetime of the BEC is about 
20 x 4.5 x 10~ 4 s = 9ms, which is close to the experiment results: the lifetime of a BEC is about 8ms. Here, by 
([9"|) , we can verify that the number of atoms in the bright soliton against the background is in the range of 4635 (at 
t = —10) and 3107 (at t = 10), which is a proper range of atoms when the soliton can be observed [10| . But, when 



4 



t 6 [-10,10], the scattering length a(t) = -0.2sech(0.197i) exp(0.02t) varies in [-0.20,-0.04] nm, which is different 
from the experiment condition: the scattering length keeps invariable when the bright soliton in BECs propagates in 
the magnetic trap (Toj |. However, up to now, they can not measure the motion of dark or bright soliton when the 
parameter varies continually, and they only measure a particular value of soliton corresponding to the fixed magnetic 
field. When we fix the magnetic field at a fixed value (the scattering length is also a fixed value) such as B — 425G 
in Ref. [1(3], the scattering length is a s = —0.21 nm, then the special value of our general solution is in accordance 
with the experimental data of Ref. [l(| • Of course, we believe that with the development of the Feshbach resonance 
technology, the experimental physicists can measure the motion of solitons in the future. Thus by modulating the 
scattering length in time via changing magnetic field near the Feshbach resonance, we may also realize the bright 
solitons in BECs. 

Fig. la Fig. lb 




IPS 51 -0.2 

-0.1-°- 15 
-0.05x 



Fig. 1: (Color online) The evolution plots of \ip u | 2 /10 4 , where A c = 4, A s = 1200, 30 = 0.4. In Fig.la: A = -0.197,7 = ki = 
-0.01; In Fig.lb: A = -0.02,7= ki = -0.001. 



When A = —0.02 (which can be derived from ujq = 2m x 7Hz and u±_ = 2ir x 710Hz) and 7 = —0.001, from Fig.lb, 
the lifetime of the BEC can reach about 200 unities of the dimensionless time corresponds to a real time of 0.1s, which 
arrives at the order of the lifetime of a BEC in today's experiments. Here, we can verify that when t is from -120 to 80, 
(i) the number of atoms is in the range of 4824 and 3234 by ([9]); (ii) the scattering length — 0.20nm < a(t) < — 0.03nm; 
(iii) From (5), we can derive £ w — — 4A^.(x — kit — xq), therefore, k\ describes the velocity of the bright 

soliton, which can be demonstrated by Fig. 1. 

It is necessary to point out that, (i) In order to give clear figures, the parameter k\ is taken to be relatively small. In 
reality, k\ is about to —4, which can be derived from the soliton's position x = —h\ exp(— At)/A; (ii) The background 
is very small with regard to the bright soliton in Fig. 1, which can be also shown by Fig. 2. Therefore the background 
may be taken as zero background approximatively; (iii) In order to keep the lifetime of the bright soliton about 8ms, 
we take 7 = —0.01 which may be the experimental value. When 7 = —0.01, by varying the value of A, it is difficult 
to extend the lifetime of the bright soliton. Thus in order to extending the lifetime of a soliton in BEC, we should 
take appropriate measures to reduce the absolute value of A and 7, as is shown in Fig. lb. 

Furthermore, we find that when cosh(£) = — 2A C cos(rj)/A s + A S /(A C cos(r])) (sinh(£) = 0), the intensity of ([5]) 
arrives at the minimum (maximum) 

,,,, ,2 , 2ri (^-4Ag)cos 2 (??) 1 sech(At) 2t 

IV>nU = A c [l- A2 _ 4 ^ cos2(7?) ]— 2— e , 

,2 _ 2 A s {Al-AAp sech(At) 2 t 
^ - [Ac + A s ~2A c cos(r,) ] ^2— e ' (8) 

This means that the bright solitons ([5|) can be only squeezed into the assumed peak matter density between the 
minimum and maximum values. Fig. 2 present the evolution plots of the maximal and minimal intensities given by 
l^ulmax/lO 4 (blue line) and iV'nlmin (y euow nne ) an d the background intensity (red line) with different parameters. 
From Fig. 2, with the time evolution, firstly the intensities increase until to the peak, then decrease to the background. 
Meanwhile, the smaller |A| and I7I, the longer the higher intensities can keep. Therefore in order to keep a bright 
soliton in BECs for longer time, we should reduce both the ratio of the axial oscillation frequency to radial oscillation 
frequency and the loss of atoms. 

To investigate the stability of the bright soliton in the expulsive parabolic and complex potential, we obtain 

L 

(\i)n\ 2 -\*Pc\ 2 )dx = N C L , (9) 
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Fig . 2a. FJ-g . 2±> 




Fig. 2: (Color online) The evolution plots of IV'n lmax/10 4 (blue line) and IV'nlmm (yellow line) and \il> c \ 2 (red line). The 
parameters in Fig. 2a-2b and Fig. 2c-2d are the same as in Fig. la and Fig. lb, respectively. 



where 

N = V a ^exp(2 7 i) (10) 

c = A a (exp(2L)-l) 

L ^l s (exp(2L) + 1) -4A c cos(77)exp(X)' 1 ' 

which is the exact number of the atoms in the bright soliton against the background described by ([5]) within [—L,L]. 
This indicates that when L takes a fixed value, for example, L = 100, then Cl — > 1, therefore the number of atoms 
in the bright soliton is determined by Nq. 
In contrast, the quantity 

L 

\^n-^ c \ 2 dx = N {C L + 4A c Mcos(r])}, (12) 



where Nq and Ci are determined by IjlOp and (fTTjl . and 

M - 1 arctanf As + 2 ^ cC0S(??) • 6Xp(L) ~ - 1 

" e arctan ^ © ex P (L) + i J ' 

9 = -4A2 C os 2 (?7), 

counts the number of atoms in both the bright soliton and background under the condition of tp(±L, t) ^ 0. Equation 
(|12|) displays that a time-periodic atomic exchange is formed between the bright soliton and the background. In the 
case of zero background, i.e., A c = 0, from (fl"2"|) the exchange of atoms depends on the sign of 7: (i) when 7 = 0, there 
will be no exchange of atoms; (ii) when 7 < 0, the exchange of atoms decreases; (iii) when 7 > 0, the exchange of 
atoms increases. As shown in Fig. 3, in the case of nonzero background and 7 < 0, a slow-fast-slow process of atomic 
exchange is performed between the bright soliton and the background, but the whole trend of the atomic exchange 
between the bright soliton and the background is decrease. In Ref. [24|, Wu et al show that the number of atoms 
continuously injected into Bose-Einstein condensate from the reservoir depends on the linear gain/loss coefficient, and 
cannot be controlled by applying the external magnetic field via Fcshbach resonance. The findings here can recover 
the same results. 

In addition, under the integration constant of J £l 4 (t)dt taken to be zero, take the following particular form at 

= (* = o,±i,±2,-,). 

tpu = fi(to) [-A c ± iA s psech{£)} exp(iA + 7*0). (13) 



This means that (fl3|) can be generated by coherently adding a bright soliton into the background. 
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Fig. 4: (Color online) The evolution plots of |i/>i 2 | 2 with A = -0.197, A c = 7,A S = 12, g = 0.4, ki = -0.01,7= -0.01. Fig. 4b 
is the contour plot. 



Inspired by two experiments [8|, [10| , we can design an experimental protocol to control the soliton in BECs near 
Feshbach resonance with the following steps: (i) Create a bright soliton in BECs with the parameters of N w 4 x 10 3 , 
u>± = 2in x 700Hz and ojq = 2n x 7 Hz, and for 7 Li. (ii) Under the safe range of parameters discussed above, 
ramp up the absolute value of the scattering length according to a(t) — — ^sech(At)exp(— 2^/t) due to Feshbach 
resonance, control the dispersion of atoms in BECs at a low level by modulating the parameter 7 about to —0.001, 
and take A to be a very small value: A = —2\ujo\/ui± = —0.02. A unity of time, At = 1 in the dimensionless 
variables, corresponds to real seconds 2/iv± = 4.5 x 10 -4 . (iii) During 200 dimensionless units of time, the absolute 
value of the atomic scattering length varies in 0.03nm < \a(t)\ < 0.20nm. This means that during the process of 
the bright soliton, the stability of soliton and the validity of ID approximation can be kept as displayed in Fig. 
lb. Therefore, the phenomena discussed in this paper should be observable within the current experimental capability. 



(II) When j3 = 0, the solution ipi is written as 



/ orzi , a 6 cosh (0 + cos(?7) + ia. sinh(Q 

^12 =U[A C + A S — r xexp(iA + 7i), (14 

cosh(£) + cos(?/) 



where AA 2 C - A 2 S > and 



1>2 
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rj = p 2 ^l 2 x — 2p 2 kx 



p\ = g (4Al-A 2 s ), 

Analysis reveals that ipi 2 is periodic with a period T — 4tt / [p 2 sech(Xt)] in the space coordinate x and aperiodic in the 
temporal variable t. Note that the period T is not a constant due to the presence of the function sech(Ai), but when 
A <C 1 and t is very small, T is very close to 4tt/p 2 . As shown in Fig. 4, when A = —0.197, A c = 7, A s = 12, g = 0.4, 
k\ = 7 = —0.01, a train of bright solitons is excited. Here the atoms in a bright soliton and in the background in 
a period [0, T] are \tpi 2 \ 2 dx ~ 3510, \A c tt\ 2 dx « 6815, respectively. Thus we can conclude that an important 
condition for exciting a train of bright solitons is that the background is strong enough. 



sech(Ai) 



B. Dynamics of a Dark Soliton in BECs 




Fig. 5: (Color online) The evolution contour plots of \ip2\ 2 , where A c — 7, A s — 12, go = 0.4, fci = —4.4, A = —0.197 and 
7 = -0.01 in Fig.5a ; A = -0.02 7 = -0.001 in Fig.5b. 

When A — * and 7 — ► 0, is reduced to the dark soliton (25| . 

iJ2 = ^y[A c + iA s tanh(^(x + Asih + ^A c ))t)]exp(iA). (15) 

Therefore, the solution (01 should be a time-dependent dark soliton, which can be shown by Fig. 5. From ([I]), we can 
obtain the intensities of the background as follows 

l^l 2 = (^ + ^)^M exp(27 i). (16) 

Therefore from (fl6|) . we guess that the solution (J4j) may describe an interesting physical process: there are "moving 
stop" , which may be realized by use of laser, at both ends of the cigar-axis direction. 
Proceeding as the case of the bright soliton, we obtain 

+00 2j4 2 

|^ 2 | 2 - \M±o°,t)\ 2 dx - -^exp(2 7 i), (17) 
-do V9o 

which describes the region of decreased density and contains a negative "number of atoms". 

As shown in Fig. 5, when the absolute of A and 7 are smaller, the dark solitons can keep a longer time and propagate 
a longer distance. Under the conditions in Fig. 5, the scattering lengths are in a range 0.05nm < a(t) < 0.20nm. 
Therefore in order to keep a dark soliton a long time in BECs, we should also reduce the values of A by adjusting the 
harmonic oscillator frequencies lo± and wq and reduce the absolute value of 7 by controlling the loss of atoms. 



III. CONCLUSIONS 



In summary, we present a direct method to obtain two families of analytical solutions for the nonlinear Schrodinger 
equation which describe the dynamics of solitons in Bose-Einstein condensates with the time-dependent interatomic 
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interaction in an expulsive parabolic and complex potential. The dynamics of a bright soliton, a train of bright 
solitons and a dark soliton are analyzed thoroughly. We can extend the lifetime of a bright soliton or a dark soliton 
in BEC by reducing the ratio of the axial oscillation frequency to radial oscillation frequency and control the loss of 
atoms. Meanwhile, our results also demonstrate that a train of bright solitons in BEC may be excited with a strong 
enough background. It is very interesting to find these new phenomena which are of special importance in the field 
of an atom laser in further experiments. 
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